num1 = conv(1.05, [0.4762, 1]);
den1 = conv(conv([0.125, 1], [0.5, 1]), [1, 1, 1]);
step(num1, den1);